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Abstract 

The Bruggeman formalism for the homogenization of particulate composite materials is used 
to predict the effective permittivity dyadic of a two-constituent composite material with one 
constituent having the ability to display the Pockels effect. Scenarios wherein the constituent 
particles are randomly oriented, oriented spheres, and oriented spheroids are numerically ex- 
plored. Thereby, homogenized composite materials (HCMs) are envisaged whose constitutive 
parameters may be continuously varied through the application of a low-frequency (dc) electric 
field. The greatest degree of control over the HCM constitutive parameters is achievable when 
the constituents comprise oriented and highly aspherical particles and have high electro-optic 
coefficients. 

Keywords: Bruggeman formalism, composite material, electro- optics, homogenization, partic- 
ulate material, Pockels effect 



1 Introduction 

In 1806, after having ascended the skies in a balloon to collect samples of air at different heights 
and after having ascertained the proportions of different gases in each sample, Jean-Baptiste 
Biot and Frangois Arago published the first known homogenization formula for the refractive 
index of a mixture of mutually inert gases as the weighted sum of their individual refractive 
indexes, the weights being in ratios of their volumetric proportions in the mixture (Biot & Arago 
1806). The Arago-Biot mixture formula heralded the science and technology of particulate 
composite materials — particularly in optics, more generally in electromagnetics, and even 
more generally in many other branches of physics. An intensive literature has developed over 
the last two centuries in optics (Neelakanta 1995; Lakhtakia 1996), and recent forays into the 
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realms of metamaterials and complex mediums (Grimmeiss et al. 2002; Weiglhofer & Lakhtakia 
2003; Mackay 2005) have reaffirmed the continued attraction of both particulate composite 
materials and homogenization formalisms. 

Post-fabrication dynamic control of the effective properties of a mixture of two constituent 
materials is a technologically important capability underlying the successful deployment of a host 
of smart materials and structures. Dynamic control can be achieved in many ways, particularly 
if controlability and sensing capablity are viewed as complementary attributes. One way is to 
infiltrate the composite material with another substance, possibly a fluid, to change, say, the 
effective optical response properties (Lakhtakia et al. 2001; Monch et al. 2006). This can be 
adequate if rapidity of change is not a critical requirement. Another way is to tune the effective 
properties by the application of pressure (Finkelmann et al. 2001; Wang et al. 2003) or change 
of temperature (Schadt & Fimfschilling 1990). Faster ways of dynamic control could involve the 
use of electric fields if one constituent material is a liquid crystal (Yu et al. 2005) or magnetic 
fields if one constituent material is magnetic (Shafarman et al. 1986). 

Our focus in this paper is the control of the effective permittivity tensor of a two-constituent 
composite material, wherein both constituent materials are classified as dielectric materials 
in the optical regime but one can display the Pockels effect (Boyd 1992). Both constituent 
materials can be distributed as ellipsoidal particles whose orientations can be either fixed or be 
completely random. 

The plan of this paper is as follows: Section 2 contains a description of the particulate compos- 
ite material of interest, as well as the key equations of the Bruggeman homogenization formalism 
(Bruggeman 1935; Weiglhofer et al. 1997) adopted to estimate the relative permittivity dyadic 
of the homogenized composite material (HCM). Section 3 presents a few numerical examples to 
show that the Pockels effect can be exploited to dynamically control the linear optical response 
properties of composite materials through a low-frequency electric field. Given the vast param- 
eter space underlying the Pockels effect, we emphasize that the examples presented are merely 
illustrative. Some brief concluding remarks are provided in Section 4. 

A note about notation: Vectors are in boldface, dyadics are double underlined. A Cartesian 
coordinate system with unit vectors ^-x,y,z is adopted. The identity dyadic is written as £, and 
the null dyadic as 0. An exp(— iuit) time-dependence is implicit with i = lo as angular 

frequency, and t as time. 

2 Theory 

Let the two constituent materials of the particulate composite material be labeled a and b. Their 
respective volumetric proportions are denoted by f a and /& = 1 — f a . They are distributed as 
ellipsoidal particles. The dyadic 

3 

= J2 a K a KS L K (1) 

K=l 
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describes the shape of particles made of material a, with ock > V-fT G [1, 3] and the three unit 
vectors ai 5 2, 3 being mutually orthogonal. The shape dyadic 
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U} b) = J2 ^b K b K (2) 

K=l 

similarly describes the shape of the particles made of material b. A low-frequency (or dc) electric 
field E dc acts on the composite material, the prediction of whose effective permittivity dyadic 
in the optical regime is of interest. 

Material a does not display the Pockels effect and, for simplicity, we take it to be isotropic 
with relative permittivity scalar in the optical regime. 

Material b has more complicated dielectric properties as it displays the Pockels effect. Its 
linear electro-optic properties are expressed through the inverse of its relative permittivity 
dyadic in the optical regime, which is written as (Boyd 1992) 



1 = i2[(l/e$+ Sj )u K u K 



K=l 



+ s 4 (u 2 u 3 + U 3 U 2 ) + S 5 (U1U3 + U3U1) + Sq (UlU 2 + U 2 Ul) , (3) 

where 
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sj=J2 r JKE d K c , J€[l,6], (4) 



K=\ 

and the unit vectors 



ui = — (u x cos 4>b + u y sin fo) cos 9b + u 2 sin 9b 
u 2 = u x sin 4>b — u y cos 4>b 
u 3 = (u x cos (f>b + u y sin cj)b) sin 9b + u z cos 6*5 



}, 6 €[O,7r], 6 €[O,27r], (5) 



are relevant to the crystallographic structure of the material. In (3) and (4), E^f = • E dc , 

K G [1, 3], are the Cartesian components of the dc electric field; ef\ 3 are the principal relative 
permittivity scalars in the optical regime; and tjk, J G [1,6] and K G [1,3], are the 18 
electro-optic coefficients in the traditional contracted or abbreviated notation for representing 
symmetric second-order tensors (Auld 1990). Correct to the first order in the components of 
the dc electric field, which is commonplace in electro-optics (Yariv & Yeh 2007), we get the 
linear approximation (Lakhtakia 2006a) 



g (b) - £ [e£> (l - e£W) u*u* 



K=l 

- 4 6) e 3 b) s 4 (u 2 u 3 + u 3 u 2 ) - ef } ef } s 5 (mu 3 + u 3 Ui) - epe^ s 6 (u x u 2 + u 2 Ui) (6) 
from (3), provided that 

{^3] 1^1} {^6]M}«1. (7) 
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This material can be isotropic, uniaxial, or biaxial, depending on the relative values of ef \ ef\ 
and e^. Furthermore, this material may belong to one of 20 crystallographic classes of point 
group symmetry, in accordance with the relative values of the electro-optic coefficients. 

Let the Bruggeman estimate of the relative permittivity dyadic of HCM be denoted by e Br . 
If the particles of material a are all identically oriented with respect to their crystallographic 
axes, and likewise the particles of material b, then e Br is determined by solving the following 
equation (Weiglhofer et al. 1997): 



fa(e (a) L-g Br ) • [l + d (a) ' (e (a) L-g Br ) 

+f b (g {b) -f) ' [l+i b) ' (g {b) -g Br 



= 0. 



In this equation, d a = d(U_^) and cf = d(uj^), where the dyadic function 



contains the unit vector 



>2vr 



sm( 



<P q =o Je q =o 



q • 



IT 1 ' e Br • IT 1 • q 



d9 n 



(8) 



(9) 



q = u 2 cos Q q + (u x cos 4> q + u y sin <f> q ) sin 9 q . 



(10) 



The Bruggeman formalism is more complicated when the relative permittivity dyadics of 
the particles of material b are randomly oriented. To begin with, let there be P > 1 distinct 
orientations. We represent the p th orientation of the relative permittivity dyadic in terms 

of the set of Euler angles {7n p }^ =1 as (Lakhtakia 1993) 

4 b) = g z ds P ) -g y {i2 p ) -g z {n P ) -cW •i; 1 (72 P ) -R- 1 ^), P = 1,2,3, ,P, 

(11) 

wherein the rotational dyadics 



| y (7) = ( u x u x + u z u z ) cos7 + (u 2 u x - u x u z ) sin 7 + u y u y 
R z {l) = {U x U x + Uj, Uy) cos 7 + (u x u y - u y u x ) sin 7 + u z u z 

Let us define 

II(7ip,72 P ,73 P ) 

Then the Bruggeman equation may be expressed in the form 



=p 



Br 



fa(e (a) L-g Br ) • [l + i a) ' ( e (a) I-| Br )] _1 + /^^S(7i P ,72 P ,73 P )=fl, 



p=l 



if all P orientations are equiprobable. In the limit P — > 00, equation (14) becomes 
fa(e (a) L-g Br ) • [L + d {a) ' (e (a) I-I Br ) 1_1 



"a~~2 / / / fl(7i,72,73) sin 72 cfyi c?72 ^73 = g. 
o7T J 73= o j 72 =o ^73=0 



-27T 
'73 



(12) 



(13) 



(14) 



(15) 
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Even more complicated orientational averages — e.g., of particulate shapes and geometric ori- 
entation of particles, in addition to crystallographic orientation — can be similarly handled. 

The HCM relative permittivity e Br can be extracted from equation (8) and equation (15) 
iteratively using standard techniques, and a Jacobi iteration technique is recommended (Michel 
2000). 



3 Numerical results and discussion 



A vast parameter space is covered by the homogenization formalism described in the previous 
section. The parameters include: the volumetric proportions and the shape dyadics of materials 
a and b; the relative permittivity scalar ; the three relative permittivity scalars and the 
upto 18 distinct electro-optic coefficients rjx of material b; the angles 9b and that describe 
the crystallographic orientation of material b with respect to the laboratory coordinate system; 
and the magnitude and direction of E dc . To provide illustrative results here, we set = 1. 
All calculations were made for two choices of material b (Cook 1996): 



I. zinc telluride, which belongs to the cubic 43m crystallographic class: 4 4 4 
8.94, r 4i = r 52 = r 63 = 4.04 x 10~ 12 m V -1 , and all other r JK = 0; and 

II. potassium niobate, which belongs to the orthorhombic mm2 crystallographic class: = 
4.72, 4 6) = 5.20, 4 b) = 5.43, r 13 = 34 x 10~ 12 m V -1 , r 23 = 6 x 10~ 12 m V -1 , r 33 = 
63.4 x 10" 12 m V -1 , r 42 = 450 x 10~ 12 m V -1 , r 51 = 120 x 10~ 12 m V -1 , and all other 
tjk = 0. 



Given the huge parameter space still left, we chose to fix /& = 0.5, the Bruggeman formalism 
then being maximally distinguished from other homogenization formalisms such as the Maxwell 
Garnett (Weiglhofer et al. 1997) and the Bragg-Pippard formalisms (Bragg & Pippard 1953; 
Sherwin & Lakhtakia 2002). Finally, we chose particles of material a and b to be spherical (i.e., 
£(a) = £(6) = ^ in g ections 3 i anc } 3 2) and spheroidal in Section 3.3. 

Two different scenarios based on the orientation of material b were investigated. The scenario 
wherein the material b particles are randomly oriented with respect to their crystallographic axes 
was considered in the study presented in Section 3.1. Particles of material b were taken to have 
the same orientation with respect to their crystallographic axes in the studies presented in 
Sections 3.2 and 3.3. 

For all scenarios, the estimated permittivity dyadic of the HCM may be compactly represented 

as 

= a Br L + p Br (u M u N + u N u M ), (16) 

wherein the unit vectors mm and u^r are aligned with the optic ray axes of the HCM (Chen 
1983; Weiglhofer h Lakhtakia 1999; Mackay & Weiglhofer 2001). For the real-symmetric rela- 
tive permittivity dyadic e Br , with three distinct (and orthonormalised) eigenvectors 61,2,3 an d 
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corresponding eigenvalues £^2,3; the scalars a Br and f3 Br are given by 



O: 



Br 



-Br 



-Br 



Br 



whereas the unit vectors um,n may be stated as 



(17) 



um 



-Br c Br \ 1 / 2 



-Br 



r Br 



-Br 



e Br ,Br 
e 3 £ 1 



ei 

1/2 



-Br ^Br\ 1 / 2 



-Br 



e3 



ei + 



-Br c Br \ 1 / 2 



e Br e Br 
e 3 e l 



e3 



(18) 



for ef r < e Br < ef\ 



In accordance with mineralogical literature (Klein & Hurlbut 1985), we define the linear 
birefringence 

(19) 



the degree of biaxiality 
and the angles 



^ = ( e r) 1/2 +( e f) 1/2 -2( e f) 1/2 



5 = cos 1 



-Br ^Br\ V 2 



On : 



cos 



e Br e Br 
e 3 e l 



COS 1 UN • U 2 



(20) 



(21) 



The linear birefringence 5 n is the difference between the largest and the smallest refractive 
indexes of the HCM; the degree of biaxiality 5^ can be either positive or negative, depending 
on the numerical value of (e^) 1 '' 2 with respect to the mean of (ef r ) 1/2 and (ef r ) 1/2 ; 25 is the 
angle between the two optic ray axes; and 6m,n are the angles between the optic ray axes and 
the Cartesian z axis. Thus, e Br can be specified by six real-valued parameters: ef r , 5 n , 5u, 5, 
and 6m,n, in a physically illuminating way 



3.1 Randomly oriented spherical electro-optic particles 

We begin by considering the scenario wherein the particles of material b are randomly oriented 
with respect to their crystallographic axes, and the particles of both materials are spherical. Ac- 
cordingly, the HCM is an isotropic dielectric medium, characterized by the relative permittivity 
dyadic g Br = e Br £. 

The Bruggeman estimate e Br , as extracted from equation (15), is plotted in Figure 1 against 
£f c , with Ef% = 0. Material b is zinc telluride for the upper graph and potassium niobate for 
the lower in this figure. The range for the magnitude of E dc in Figure 1 — and for all subsequent 
figures — was chosen in order to comply with (7). In the case where material b is zinc telluride, 
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e varies only slightly as changes, and e is insensitive to the sign of e . A greater 
degree of sensitivity to E$ c is observed for the HCM which arises when material b is potassium 
niobate; in this case the HCM's relative permittivity is also sensitive to the sign of E^ c , thereby 
underscoring the significance of crystallographic class of the electro-optic constituent material 
even when averaging over crystallographic orientation is physically valid. 

3.2 Identically oriented spherical electro-optic particles 

The scenario wherein all particles of material b are taken to be identically oriented, with particles 
of both constituent materials being spherical, is now considered. Let us begin with the case 
where material b is zinc telluride. This is an isotropic material when E dc = 0, but the anisotropy 
underlying the Pockels effect becomes evident on the application of the low-frequency electric 
field (Lakhtakia 2006b). 

The HCM parameters, as extracted from equation (8), are plotted in Figure 2 as functions 
of Ef c 2 with E^ c = 0. The crystallographic orientation angles 6b = (fib = 0. As expected, in this 
figure 5 n = 5u = (i.e., the HCM is isotropic) when E( c = E 2 C = 0. The HCM constitutive 
parameters e 2 r , S n , 5u, 5, and 9m,n are all insensitive to the signs of E( c and Ef. The HCM 
is negatively biaxial in general (because Su < 0), although the biaxiality is small. The linear 
birefringence is not sensitive to the signs of Ef c 2 ; it increases considerably as |E dc | is increased, 
as a glance at data on minerals readily confirms (Griblle & Hall 1992). The two optic ray axes 
remain almost mutually orthogonal, as indicated by 5 ~ 45°, as \Ef c 2 \ are changed. 

The influence of the orientation angle 9b for zinc telluride is explored in Figure 3. Here, the 
optic ray axis angles 9m,n are plotted as functions of Ef c 2 with E$ c = 0, for 9b £ {45°, 90°} 
with 4>b = 0. The orientations of both optic ray axes continuously vary with increasing Ef c 2 
in a manner which continuously varies as 9b increases. The polar angle of the optic ray axis 
aligned with um, namely 9m, is slightly sensitive to E 2 C but insensitive to Ef c . In contrast, 
9n is acutely sensitive to both Ef c and E 2 C . Furthermore, 9n is sensitive to the sign of E 2 C 
but not the sign of Ef c . The HCM parameters e 2 r , S n , Sbi, 6 — which are not presented in 
Figure 3 — are insensitive to increasing 9b; the plots for these quantities are not noticeably 
different from the corresponding plots presented in Figure 2. 

Let us turn now to the case where material b is potassium niobate. This material is anisotropic 
(orthorhombic and negatively biaxial) even when the Pockels effect is not invoked, and it has 
much higher electro-optic coefficients than zinc telluride — hence, it can be expected to lead a 
different palette of HCM properties. 

The HCM parameters are plotted in Figure 4 as functions of E 2 \ with Ef c = 0. As in Figure 2, 
the crystallographic orientation angles of material b are taken as 9b = (fib = 0. Whereas the 
parameters e 2 r , 5 n , Sbi, and 5 are not particularly sensitive to E 2 C , they do vary significantly as 
E^ c varies. Most notably, the HCM can be made either negatively biaxial (Sbi < 0) or positively 
biaxial (Sbi > 0)- The two optic axes of the HCM need not be mutually orthogonal, with the 
included angle 26 between them as low as 40°. The polar angles 9m,n are sensitive to both E 2 C 
and E^ c . We note that the sign of Ef 2 does not influence any of the HCM parameters, but the 
sign of E 2 C does influence the polar angles 9m,n- 

The influence of the orientation angle 9b is explored in Figure 5. The constitutive parameters 
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of material b are the same as in Figure 4 but with 9b £ {45°, 90°}. As is the case in Figure 3, the 
graphs in Figure 5 show that the dependencies of the polar angles 9m, N upon the components 
of E dc are acutely sensitive to 9b- The HCM parameters e^ r , 5 n , 6^, and 5 — which are not 
presented in Figure 5 — are insensitive to increasing the plots for these quantities are not 
noticeably different to the corresponding plots presented in Figure 4. 

A comparison of Figures 2 and 3 with Figures 4 and 5 shows that the application of E dc is 
more effective when material b is potassium niobate rather than zinc telluride. A dc electric 
field that is two orders smaller in magnitude is required for changing the HCM properties with 
potassium niobate than with zinc telluride, and this observation is reaffirmed by comparing the 
upper and lower graphs in Figure 1. To a great extent, this is due to the larger electro-optic 
coefficients of potassium niobate; however, we cannot rule out some effect of the crystallographic 
structure of material b, which we plan to explore in the near future. 

Electrical control appears to require dc electric fields of high magnitude. However, the needed 
dc voltages can be comparable with the half-wave voltages of electro-optic materials (Yariv &; 
Yeh 2007) . We must also note that the required magnitudes of E rfc are much smaller than the 
characteristic atomic electric field strength (Boyd 1992). The possibility of electric breakdown 
exists, but it would significantly depend on the time that the dc voltage would be switched 
on for. Finally, the non-electro-optic constituent material may have to be a polymer that can 
withstand high dc electric fields. 

3.3 Identically oriented spheroidal electro-optic particles 

We close by considering the scenario wherein the effect of the Pockels effect is going to be highly 
noticeable in the HCM — that is, when the particles of material b are highly aspherical and the 
crystallographic orientation as well as the geometric orientation of these particles are aligned 
with E rfc . We chose potassium niobate — which is more sensitive to the application of E rfc than 
zinc telluride — for our illustrative results. 

In Figure 6, the HCM parameters ef r , 5 n , 5^, 6, and 9m are plotted against E$ c , with Ef^ = 
0. Both constituent materials are distributed as identical spheroids with shape parameters 
ai 5 2 = Pi,2 = 1 and 03 = 03 £ {3,6,9}; furthermore, 9b = <fib = 0. As 9n = 180° — 9m for 
this scenario, 9n is not plotted. All the presented HCM parameters vary considerably as Ef 2 
increases; furthermore, all are sensitive to the sign of E^ c . 

We note that the degree of biaxiality and the linear birefringence increase as a.3 = increases. 
This is a significant conclusion because perovskites (such as potassium niobate) are nowadays 
being deposited as oriented nanopillars (Gruverman & Kholkin 2006). 

4 Concluding remarks 

The homogenization of particulate composite materials with constituent materials that can 
exhibit the Pockels effect gives rise to HCMs whose effective constitutive parameters may be 
continuously varied through the application of a low-frequency (dc) electric field. Observable 
effects can be achieved even when the constituent particles are randomly oriented. Greater 
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control over the HCM constitutive parameters may be achieved by orienting the constituent 
particles. By homogenizing constituent materials which comprise oriented elongated particles 
rather than oriented spherical particles, the degree of electrical control over the HCM consti- 
tutive parameters is further increased. The vast panoply of complex materials currently being 
investigated (Grimmeiss et al. 2002; Weiglhofer & Lakhtakia 2003; Mackay 2005; Mackay & 
Lakhtakia 2006) underscores the importance of electrically controlled composite materials for a 
host of applications for telecommunications, sensing, and actuation. 
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Figure 1: The estimated relative permittivity scalar e Br of the HCM plotted against E^ c (in 
V m _1 xlO 8 ) for Ef c 2 = and f a = 0.5. Material b is zinc telluride for the upper graph and 
potassium niobate for the lower graph. The particles of material b are randomly oriented with 
respect to their crystallographic axes, and both types of particles are spherical. 
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Figure 2: The HCM parameters , S n , 5u, 5 (in degree) and 9m,n ( m degree) plotted against 
Ef°2 (in V m _1 xlO 9 ). The crystallographic orientation angles of material b are 9^ = 4>b = 0> 
and £f c = 0. Material b is zinc telluride, and the particles of both constituent materials are 
spherical. 
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Figure 3: The optic ray axis angles 9m,n (hi degree) of the HCM plotted against Ef c 2 (in V 
m _1 xlO 9 ). The constitutive parameters of material b are same as in Figure 2 except that 
9b S {45°, 90°}. (The corresponding plots of ef r , S n , 5u, and 5 are not noticeably different to 
those presented in Figure 2.) 
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Figure 4: The HCM parameters ef r , 8 n , 8u, 8 (in degree) and Bm,n (in degree) plotted against 
£"23 (in V m _1 xlO 7 ). The crystallographic orientation angles of material b are 6b = 4>b = 0; 
and Ef c = 0. Material b is potassium niobate, and the particles of both constituent materials 
are spherical. 
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Figure 5: The optic ray axis angles 6m,n (in degree) plotted against £"23 (in V m _1 xlO 7 ). 
The constitutive parameters of material b are same as in Figure 4 except that 9b £ {45°, 90°}. 
(The corresponding plots of 63 , <5 n , &Ui and <5 are not noticeably different to those presented in 
Figure 4.) 



16 



12 3 



Figure 6: The HCM scalar parameters ef r , S n , and 5u, 5 (in degree), and 9 m (in degree) plotted 
against E^ c (in V m _1 xlO 8 ). The crystallographic orientation angles 9b = <pb = 0, and Ef c 2 = 0. 
Material b is potassium niobate. The constituent materials are distributed as spheroids with 
shape parameters: a\ = ct2 = (3\ = fa = 1; and a% = = 3 (solid curves), 03 = (3^ = 6 
(dashed curves), and «3 = 0s = 9 (broken dashed curves). The optic ray angle 9n = 180° — 9 m- 
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